arXiv:1509.01235vl [gr-qc] 3 Sep 2015 


Macroscopic Einstein Equations for a Cosmological Model] Macroscopic Ein¬ 
stein Equations for a Cosmological Model with A - Term 

Abstract 

Through averaging the Einstein equations over transverse gravitational 
perturbations it is obtained a closed system of two ordinary differential 
equations describing macroscopic cosmological evolution of the isotropic 
space-flat Universe filled with gravitational radiation. It is found an 
asymptotic solution of evolution equation for gravitational perturbation 
amplitude. Making the substitution of this solution into Einstein equation 
averaged over gravitational perturbations, the single evolution non-linear 
ordinary differential second-order equation relative to macroscopic scale 
factor is obtained. It is also found a solution of evolution equation for scale 
factor in WKB-approximation which analytically describes the process of 
transformation from ultrarelativistic regime of cosmological extension to 
inflationary one. 

Keywords: macroscopic Einstein equations, transverse gravitational per¬ 
turbations, evolutions equations, inflation 

1 Gravitational Perturbations of Isotropic Uni¬ 
verse 

Gravitational perturbations of homogenous isotropic Universe were first intro¬ 
duced and classified in pioneer work by E.M. Lifshitz, 1946 [T] (see e.g. the 
well-known L.D. Landau and E.M.Lifshitz monograph m)- In classical works of 
R.A. Isaakson, 1968 [31H] there was developed the approach to build macroscopic 
gravitation theory based on averaging of microscopic (short-wave) fluctuation 
of metrics. In Author’s works 1985-1991 (see e.g. |S], 2007) it was being devel¬ 
oped the statistical theory of gravitational interaction based on combination of 
ideas expressed in R.A. Isaakson’s and N.N. Bogolyubov’s works. In particular, 
these ideas were realized in [^ for dynamical derivation of a kinetic equation 
for photons against a background of locally fluctuating but macroscopically ho¬ 
mogenous isotropic Universe. In the article we realize these ideas for a case 
of macroscopically flat Universe described by Einstein equations with a cosmo¬ 
logical term assuming that single matter type in this Universe is presented by 
gravitational waves i.e. transverse traceless gravitational perturbations. As is 
well known, in absence of these perturbations the Universe is described by de 
Sitter solution or, if using synchronous system of coordinates - by inflationary 
solution. Therefore not only process of derivation of closed evolution equations 
describing macroscopic Universe is of interest for us but also the research of 
such dynamically reasoned macroscopical model’s behavior. 
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1.1 Averaging of Einstein Equations 

According to general approach to derivation procedure of macroscopic equations 
[3 let us represent space-time metrics in form: 

9ik — 9ik A ^9ik — ^9ik (1) 

where e <C 1 is a macroscopic metric of space-time obtained using certain 
averaging: 

gfk = Wk^ ( 2 ) 

so that: 

6gik = 0 . ( 3 ) 

Assuming averaging operation to be independent on coordinates, let us also 
require derivatives of all arbitrary perturbations averages to be equal to zero: 

djSgik = 0 ; djiSgik = 0 . ( 4 ) 

Let us consider Einstein equations with cosmological tern|3: 


Gik ^gik — 0 , 


(5) 


where Gik = Rik — ^/‘^Rgik is Einstein tensor. 

Let us write Einstein equations with cosmological term in the second ap¬ 
proximation over gravitational perturbations: 

Gik + = ^9ik, ( 6 ) 


where Gf^} = G^ugfk , G^^} = Un[Gikidg,k)] eG[l\ G^f = Uii[Gik{S'^gik)] 

and average these equations taking into account relations (O, Q: Thus 
we get macroscopic second-order Einstein equations over perturbations of grav¬ 
itational field: _ 


/^(O) _ _^(2) 

^ik — ^ik 


XJO) 
^9ik ’ 


(7) 


according to which second-order corrections can be considered as gravitational 
perturbations’ energy-momentum tensoid: 


T^k = -^dS. ( 8 ) 

^In the article Ricci tensor is defined through contraction of curvature tensor over first and 

third indices, metrics signature is (-h). 

^This statement is one of the main statements of R.A. Isaakson’s theory mm 
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1.2 Transverse Perturbations of Space-Flat Isotropic Uni¬ 
verse 


Let us write a metrics with monochromatic transverse gravitational perturba¬ 
tions of space-flat isotropic Universe in form (see e.g. my- 


dsQ = a^{f]){dr]'^ — dx^ — dy^ — dz'^)-, 

(9) 

ds^ = dsQ -f a? {ifjhapdx^dx^ 

( 10 ) 

hai3 = eai3S{r])e"’^'^, 

(11) 

where S{r]) is an amplitude of gravitational waves. Thus, it is: 


d'lk = a^(^)Diag(-l, -1, -1, +1); 

( 12 ) 

^94a — O 7 ^ 

(13) 

Then: 


hp = 

(14) 

h = /i“ = go^hap = + ^22 + hss), 

(15) 

herewith for transverse perturbations it is: 


hpUa = 0 ; 

(16) 

h = 0. 

(17) 

As a result of (11711 in approximation linear over /i it is: 



(18) 


by all directions of wave vector 


1.3 Averaging Perturbations over Wave Vector Direction 

Thus, we expand Einstein tensor into series in smallness of gravitational waves 
S{r]). Herewith using isotropy of unperturbed metrics it is convenient to intro¬ 
duce local system of coordinates where : 

n = n(0,0,l); s = (l,0,0), (19) 

where s is a unit vector of transverse perpetrations’ polarization. In this system 
of coordinates it is 


hi2 = 0; hn = -h22 = ^(r/)e“^ (20) 

^It is also possible to carry out averaging over all lengths of wave vectors but this operation 
does not provide additional information. See 013 about averaging procedure and getting 
Einstein macroscopic equations . 
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In arbitrary Cartesian coordinate system of 3-dimensional Euclid space -E 3 po¬ 
larization tensor Cik of formula (I15|) has a form: 

_ o I . , 

ea/3 — 9 ^a0i V^f) 

§2 = 1; sn = 0, xv^=r?. (22) 

It is easily checked that calibration condition (ITSl) is automatically fulfilled in 
such a case. 

Let us note that against a background of isotropic space, operation of av¬ 
eraging over directions is reduced to calculating an integral over 2 -dimensional 
sphere of radius n: 

()i(n, r) = ^ J (23) 

Thus, according to ©, ® we have: 

0, 0,... (^4) 


1.4 Zero Approximation 

Expanding Einstein tensor in metrics perturbation, we get known expressions 
in zero approximation: 


^(0) _ ^(0) _ ^(0) _ 

1^11 _ (^22 — ^33 ~ ^ ^ ^2 ’ 

/2 


^(0) _ _ 

^44 — ^ o • 


, a 


(25) 

(26) 


Thus, in zero over gravitational perturbations approximation, we could get stan¬ 
dard Einstein equation with A - term 


and its inflationary solution: 



(27) 


(28) 


1.5 Equation for Gravitational Wave Amplitude 

In approximation linear over S we get single non-trivial components: 


G-W = -G' 22 ’ =SG= ^e*"^ 


Ai) 


1 


S" + 2-S' + s(n^ + 2^ - 4— 
a \ a 


(29) 


Covariantly generalizing the result in i? 3 , let us write: 

— 2SaSp)5G. 


(30) 
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Substituting expression (1^ in Einstein equations ([5]), we get equation for 
an amplitude of transverse perturbations: 

S" + 2-5" + s(n'^+ 2^-4— + = 0. (31) 

a \ a j 

This equation with an account of relations (I25|) can be written in simplified 
form: 

5" + ^5' + (n^ - 2G^°^ + \a^)S = 0. (32) 

In particular, substituting zero approximation inflationary solution of Einstein 
equations ([28|l here, equation (1^ can be reduced to the next one: 

+ (33) 

which has its solution: 

5 = Giry3/2j^(„^) + (34) 

M=^\/3 + 4/A. (35) 

Particularly, near cosmological singularity of zero approximation ry —> — oo, 
hence |nry| —>• oo equation (|32|) is reduced to the next oscillation equation: 

5" + n^S = 0, (36) 

having its solution ordinary WKB-solutions: 

5 = Gie™’’ + ^ 26 "“’^, (37) 

which can be obtained also from the exact solution of (IMl) in this limit. 


1.6 Second Approximation 

Calculating Einstein tensor of second approximation we obtain its non-trivial 
components: 


Gii = G^22 = ’ 

= e^™^ f-S^n^ + 55" + -5'^ + 2-55'^, 
V4 4a/ 


= -e^ 


2inz 


-S^n^ + ^5'^ + 2-55' ). 
4 4a 


Carrying out covariant generalization in Ea, we can write: 


(38) 

(39) 

(40) 


(41) 
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where 


(42) 

(43) 


u = + SS" + + 2-SS'] 

4 4 a 

y = S'2„2_ 1^/2^ 

Averaging (1411) over directions of perturbation propagation with an account of 
obvious equality 

naUp = (44) 

we get for averaged components the next expression: 

= 87r(£: + V)u^uk - 8^TVg^k , (45) 


where u* is a timelike vector of observer velocity and £ and V are energy density 
and pressure of transverse gravitational perturbations: 


£ = 




- 2-SS' 

4 4a 


V = - 


J_(lls^n^ + -S’^ + 2SS'- 

Stto? V12 12 a 


SS" 


(46) 

(47) 


In particular, for WKB-solution (1571) these formulas lead to effective ultra- 
relativistic state equation: 


V 


£ 
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1 

167ra^ 


5^ 



(48) 

(49) 


2 Macroscopic Einstein Equations of Second Or¬ 
der of Gravitational Perturbations for Isotropic 
Space-Flat Universe 

2.1 Evolution Equations 

Combining obtained results within the frame of equations ([5]) and ©, we ob¬ 
tain a self-consistent system of ordinary second-order differential equations de¬ 
scribing cosmological evolution of flat macroscopic Universe with an account of 
transverse gravitational perturbations: 

S" + 2-S’ + s(n^ + 2^^ - A— + \aA = Q- (50) 

a \ a ) 

3^^ =-h ^5'^ - 2-S'S"-k Aa^ (51) 

4 4 a 
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In such a case equation (ISO)) describes cosmological evolution of scalar amplitude 
S{r]) of gravitational perturbations and equation (ISTj) describes cosmological 
evolution of scale factor 0 ( 77 ). In this case evolution equation for gravitational 
perturbations is linear differential equation of second order relative to amplitude 
of these perturbations. 

Let us notice that if there is a perturbation spectrum instead of monochro¬ 
matic perturbations (El): 

J e„Mn) + 5:(r7)e—’] d^n, (52) 

then in evolution equation (l50l) it is necessary to make a substitution S ^ Sn 
and in evolution equation for scale factor m it is necessary to use expressions 
for averages: 


—>■ n'^S^ = 

S'2 = 


1 


v^. 

I 


V^, 


S'S S'S = 


1 


1 , 




SnS^n^dri] 

(53) 

r*oo 

S'nS*n'dn-, 

(54) 

+ Sr,S:,')dn. 

(55) 


2.2 Asymptotic Solutio of Evolution Equation for Pertur¬ 
bations 

Let us notice that using scale transformation of amplitude S 



a 

in equation (l50l) it is possible to get rid of first derivative: 

(/>" + Q{n,r))(j) = 0, 


where 

Q{n, r]) = n'^ + 2 ^ _ 5 ^ _|_ 
a 

Let us consider asymptotics of solutions of this equation assuming 


(56) 


(57) 

(58) 


n » I; Q{n, ij) > 0. (59) 

Using theorem [5] about asymptotic solution of equation (EZl), let us write its 
asymptotic independent solutions: 

^ exip^^±i y = Q(n, (60) 

^0 
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( 61 ) 


and asymptotic values of derivatives (taking into account equations 

T) 

(j)'{r]) ±iQ{n,r]y^^expi^±i j yjQ{n,r]')dr]''^ = zbi(5(n, 

r)o 

77 

^ -Qin,r]f^^expi^±i J \/Qin,r]')drj'j = (62) 

VO 

Further, for us it is important the fact that terms which are, practically, 
quadratic over amplitudes and included in the right part of evolution equation 
(EOl) are calculated as a square of complex variable’s modulus i.e. = SS*, 
where S and S* are two independent asymptotics O- Therefore it is: 


S^^SS* = ^r-, S'^ ^ S'S*' = 
- (# + (j> (!>*)— +4>(j) — = 

1 


SS' ^ + 55*') = + #*') - 

Calculating, we find: 


52 ~ 


02 

Oq 


a‘^{v)\/Q(n^ 


55'~- 


o2 

Oq 


a^irhyjQin.r]) a ’ 


gl2 ^ g2 \/Q('»-, Vi) 


-^0 Y 


a‘^{'n)\/Q{n,ri) 


(63) 


(64) 


where 5o is a constant. 

Substituting these expressions in equation of scale factor evolution (IFTl) . we 
finally get macroscopic evolution equation in asymptotic approximation: 


3^-A=4fI 


a V4 y(5(n,ry) 4 


+ + 7 


9 a'- 


4 \/Q{n,ri) 


(65) 


hus, we obtained a closed ordinary highly non-linear differential second-order 
equatioifl, which provides an asymptotically exact description of macroscopic 
cosmological evolution of the Universe filled with weak transverse gravitational 
perturbations. Basically, this equation can be investigated by methods of quali¬ 
tative theory of differential equations. We intend to revisit this problem in next 
publications. Now let us consider WKB-approximation of equation (1651) . 


2.3 WKB-Solution of Evolution Equation for Scale Factor 

Let us consider now the following WKB-approximation of evolution equations: 

a' 9 a" 

n » —; n — => nry 1. (66) 

a a 

^let US remind that function Q{n,rf) depends on a,a',a'' 




















In this approximation 


Q{n,r]) Ri n, 

(67) 

equation (1^51) takes fairly simple form: 


a'2 2Sin 

a4 a* 

(68) 


Solution of this equation can be written in elementary functions, proceeding to 
physical time t, so that adrj = dt. Carrying out elementary integration, we find: 


a{t) = 



(69) 


Solution of such form was obtained by the Author (see e.g. [9]). This solution 
describes smooth transition from ultrarelativistic stage of Universe expansion 
to inflationary stage a.t t > tc = y^3/4A. 

Actually, at —>■ 0 we get from ([69|) : 



and at t —>■ oo we get from (IMl): 


a{t) 




(70) 


(71) 


3 Discussion of Results 

Summarizing the results of the article, let us indicate the most important ones: 

1. Equations of second order over perturbation amplitude are obtained from 
Einstein equations with A - term using expansion of metrics in small trans¬ 
verse perturbations relative to background Friedmann solution. 

2. Eventually, a closed mathematical model describing cosmological evolu¬ 
tion of the Universe filled with gravitational radiation, is obtained. This 
model consists of two ordinary second-order differential equations which 
we call evolution equations for the sake of brevity. 

3. First equation for evolution of amplitude of gravitational perturbations’ 
monochromatic mode is a linear homogenous differential equation. 

4. Second evolution equation describes evolution of macroscopic scale factor 
of the Friedmann Universe. This equation is highly non-linear and is de¬ 
fined trough solutions of evolution equation for amplitude of gravitational 
perturbations. 
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5. An asymptotic solution of evolution equation for amplitude of gravita¬ 
tional perturbations is found. 

6. Using obtained solution there have been calculated the macroscopic aver¬ 
ages of perturbation amplitude’s square and its derivatives. 

7. Eventually, closed highly non-linear ordinary differential equation describ¬ 
ing cosmological evolution of scale factor is obtained. Energy spectrum of 
gravitational perturbations and cosmological constant are parameters of 
this equation. 

8. It is found the WKB-solution of this evolution equation analytically de¬ 
scribing transformation from ultrarelativstic stage of the Universe expan¬ 
sion to inflationary one. 

Thus, as was noticed in |10j , an account of gravitational perturbations in energy 
balance of the Universe at its early stages really discards primordial inflationary 
expansion stage and puts it on the second place after ultrarelativstic stage. 

Let us notice that in recent work Chiu Man Ho and Stephen D.H. m, de¬ 
voted to quantum instability of the de Sitter Universe due to particles generation 
it is obtained similar conclusion. 

Further we will be purposing first, on investigation of solution of evolu¬ 
tion equations by numerical methods and second, on construction of equivalent 
mathematical model accounting other types of gravitational perturbations and 
physical fields. 
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